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Haydeeite: a spin-1/2 kagome ferromagnet 
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The mineral haydeeite, <a-MgCu 3 (OD) 6 Cl 2 , is a S = 1/2 kagome ferromagnet that displays long- 
range magnetic order below Tc = 4.2 K with a strongly reduced moment. Our inelastic neutron 
scattering data show clear spin-wave excitations that are well described by a Heisenberg Hamiltonian 
with ferromagnetic nearest-neighbor exchange Ji = — 38 K and antiferromagnetic exchange Jd — 

+ 11 K across the hexagons of the kagome lattice. These values place haydeeite in the ferromagnetic 
part of the phase diagram, close to the non-coplanar twelve-sublattice cuboc2 antiferromagnetic 
order. Diffuse dynamic short-range ferromagnetic correlations observed above Tc persist well into 
the ferromagnetically ordered phase with a behavior distinct from critical scattering. 

PACS numbers: 75.10.Kt, 75.10.Jm, 75.40.-s, 75.40.Gb, 78.70.Nx 


Quantum S = 1/2 spins localized on the two- 

dimensional kagome lattice of vertex-sharing triangles are 
highly frustrated magnetic model systems that are pre¬ 
dicted to form quantum spin liquids (QSLs) with exotic 
ground states and excitations pQ. While initial work fo¬ 
cused on nearest-neighbor Heisenberg kagome antiferro- 
magnets EHZ], such as the mineral herbertsmithite [8], 
it has been realized that further neighbor interactions 
or the antisymmetric Dzyaloshinski-Moriya (DM) inter¬ 
action, which is intrinsically allowed in the kagome lat¬ 
tice, may lead to even more fascinating ground states 
[MS. Even for the case of classical spins, the kagome 
lattice with interactions beyond nearest neighbors (NN) 
[see Fig. [lja)] shows a very rich magnetic phase diagram 
with complex non-coplanar structures such as chiral and 
multi-k magnetic ordering [IMS]. This is illustrated 
by the uniformly filled areas in Fig. [ljb) for the case of 
ferromagnetic NN interaction Ji < 0. Very recently, it 
was shown that for quantum spins, QSL phases “invade” 
parts of the classical phase diagram [151118] , as shown by 
the open circles in Fig. [ljb). 

An experimental realization of such a system that ex¬ 
emplifies the importance of quantum effects is the min¬ 
eral kapellasite [T9H23] , where high-temperature fits to 
the magnetic susceptibility [22 place it in the area of 
the phase diagram that corresponds to the non-coplanar 
twelve-sublattice magnetically ordered cuboc2 phase [see 
(blue) solid square in Fig. [ljb)]. However, kapellasite 
does not actually order magnetically, but instead forms a 
QSL with gapless spinon-like excitations and short-range 
correlations, which are reminiscent of the cuboc2 phase 
[2T] . Following recent theoretical work m , it can now 
be understood that kapellasite falls in the region of the 
phase diagram where quantum fluctuations favor a chi¬ 
ral gapless spin liquid state labelled CSL-A, which has 


similar short-range correlations to kapellasite. 

In this Rapid Communication, we explore experimen¬ 
tally other parts of the phase diagram shown in Fig.CTb). 
By replacing the Zn 2+ ion [green open circle in Fig. Ilia)] 
of kapellasite by Mg 2+ , we expect the Jg exchange inte¬ 
gral to change the most. Such a substitution leads to the 
mineral haydeeite, which orders ferromagnetically m • 
Our inelastic neutron scattering data allow us to charac¬ 
terize the ground state and determine the exchange inte¬ 
grals via spin-wave theory. We find that it is the decrease 
of J d /1 J\ | that places haydeeite in the ferromagnetic part 
of the phase diagram rather than in the cuboc2 antifer- 
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FIG. 1. (Color online) (a) Kagome lattice of haydeeite with 
Cu 2+ 5 = 1/2 spins (blue solid circles), diamagnetic Mg 2+ 
ions (green open circles), and main exchange interactions (red 
dashed lines), (b) Phase diagram of the kagome lattice for a 
ferromagnetic nearest-neighbor interaction, J\ < 0. The filled 
(colored) areas correspond to magnetically ordered phases and 
the open circles show regions where quantum spin liquid states 
are predicted [18]. The (blue) solid square represents kapell¬ 
asite 21 , [22] while the (red) solid circle corresponds to hay¬ 
deeite. 
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FIG. 2. (Color online) Magnetic Bragg scattering from hay- 
deeite, <a-MgCu 3 (OD) 6 Cl 2 , obtained by integration over the 
elastic peak (|F?| < 0.02 meV) of data from IN5 taken at 
T — 1.5 K after subtraction of 8 K data. The (red) dots are 
experimental data with one-sigma error bars, the solid (black) 
line going through the data points is the calculated diffraction 
pattern, and the (blue) line at the bottom is the difference. 
The inset shows the inverse susceptibility with a linear fit. 


romagnetic state. 

Haydeeite, <a-MgCii 3 (OH) 6 Cl 2 , crystallizes in the trig¬ 
onal P3ml space group (#164) with lattice parameters 
a = 6.2885(1) and c = 5.7271(1) A. It is isostructural 
with kapellasite mm and also has 5=1/2 Cu 2+ spins 
forming a highly two-dimensional kagome lattice of equi¬ 
lateral corner-sharing triangles. The kagome lattice is 
sketched in Fig. [lja) with the main exchange interac¬ 
tions Ji, J 2 , and Jd. The diamagnetic Mg 2+ ions sit at 
the centers of the hexagons. Note that haydeeite should 
not be confused with its polymorph ‘Mg-herbertsmithite’ 
7 -MgCu 3 (OH) 6 Cl 2 , which has the same chemical formula 
but a different crystallographic structure and different 
magnetic properties [24] . 

A 1.75 g powder sample of fully deuterated haydeeite 
was synthesized using a scaled-up variant of the estab¬ 
lished method dang. Powder diffraction using x-rays 
and neutrons confirms the crystallographic structure [20] 
and, as expected, a higher degree of site order than in 
kapellasite. In haydeeite, there are no Cu ions on the 
Mg site and about 16% Mg ions reside on the Cu sites, 
resulting in a slightly diluted kagome lattice. 

Measurements of the magnetic susceptibility x(T) at a 
field of H = 5 T and the magnetization M(H) at T = 2 K 
show that deuterated and hydrogenated haydeeite have 
very similar magnetic properties: Ferromagnetic ordering 
sets in at Tc = 4.2 K, the saturated moment is 0.83 fiB 
per mole of Cu, the Curie-Weiss temperature extracted 
from a fit of the linear region 150 < T < 300 K of the 
inverse susceptibility is #cw = 28 ± 3 K, and a drop in 


\T at low temperatures suggests antiferromagnetic com¬ 
ponents. 

Inelastic neutron scattering measurements on the fully 
deuterated <a-MgCu 3 (OD) 6 Cl 2 powder sample were per¬ 
formed on the time-of-flight (TOF) spectrometers IN5 
and IN4 at the Institut Laue-Langevin, using neutrons 
with incident energies Ei of 3.55 and 16.6 meV. These 
measurements show both acoustic and optic spin-wave 
branches in the ferromagnetically ordered phase as well 
as weak magnetic Bragg peaks that are apparent in the 
difference data (Fig. |2|. The high sensitivity, low back¬ 
ground, and energy selectivity of the inelastic TOF tech¬ 
nique allows us to conclude that these Bragg peaks, which 
were not observed by standard (energy-integrating) neu¬ 
tron diffraction techniques, probably due to the overlap¬ 
ping strong inelastic diffuse scattering discussed below, 
correspond to a k = 0 propagation vector. Combined 
with the magnetic susceptibility measurements, this con¬ 
firms the ferromagnetic ordering. 

However, the inelastic TOF technique is not well suited 
to actually determine the magnetic structure. We have 
compared on a qualitative level the observed magnetic 
intensity to that calculated for the three irreducible rep¬ 
resentations (IRs) [2E]. The best description is obtained 
with a collinear magnetic structure where the moments 
are in the kagome planes, which corresponds to the two- 
dimensional T 5 IR in Kovalev’s notation [27j, see Fig. 
[ 2 ] A rough estimate of the magnetic moment gives 
a strongly reduced magnitude of at most 0.2 (Ib at 
T — 1.5 K, which suggests strong quantum fluctuations. 

The observed inelastic neutron scattering intensity 
is proportional to the dynamic magnetic susceptibility 
x"(Q,E), shown in Fig. [3] for <a-MgCu 3 (OD) 6 Cl 2 before 
correction for the magnetic form factor. Acoustic spin 
waves are clearly observed in the low-energy data of the 
ferromagnetically ordered phase shown in Fig. [3ja), and 
are even more apparent after subtraction of the param¬ 
agnetic high-T data of Fig.|3jb), see Fig.[3|c). Optic spin 
waves are observed at higher energies, see Fig. [3je). 

The magnetic scattering can be modeled using stan¬ 
dard linear spin-wave theory for an isotropic Heisenberg 
Hamiltonian with further-neighbor interactions. Since 
haydeeite has three magnetic ions per primitive mag¬ 
netic unit cell and a propagation vector km 0 , there 
are three spin-wave modes, whose energies and intensities 
are obtained by diagonalizing a real 3x3 matrix. Stan¬ 
dard Monte Carlo techniques are then used to obtain the 
powder-averaged dynamic susceptibility x"(Q,^). The 
best agreement with the data is found for J\ = —38, 
J 2 = 0, and Jd = +11 K, where negative exchange inte¬ 
grals correspond to ferromagnetic interactions and each 
bond is counted once in the Hamiltonian. The corre¬ 
sponding calculated x"(Q,E) is shown in Figs.[3^d) and 
if)- The quality of the agreement is clearly seen in 
Q cuts at constant energy, see Fig. [4] The calculated 
spectra are very sensitive to the value of J 2 , and from 
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FIG. 3. (Color online) Dynamic susceptibility x\Q^E) on a 
linear intensity scale as a function of wave vector Q and energy 
E of fully deuterated haydeeite. Low-energy data measured 
with Ei = 3.55 meV are shown in the panels: (a) ferromag¬ 
netic state at T = 1.5 K, (b) paramagnetic state at T = 8 K, 
and (c) magnetic part of the scattering after subtraction of 
the T — 8 K data. Panel (e) shows the magnetic scattering 
at higher energies (Ei = 16.6 meV) in the ferromagnetic state 
at T = 2.2 K after subtraction of the high-temperature data 
at T — 60 K. Spin-wave calculations at low and high energies 
are shown in panels (d) and (f), respectively. 


our spin-wave analysis we infer that |J 2 |/|Ji| << 0.1. 
The Curie-Weiss temperatures calculated from these ex¬ 
changes is #cw = 32 K, in close agreement with the 28 K 
estimated from the inverse susceptibility. 

Density functional theory (DFT) has predicted a rel¬ 
atively strong antiferromagnetic nearest neighbor inter¬ 
action J\ in both haydeeite and kapellasite [28], in con¬ 
trast to the ferromagnetic J\ observed experimentally. 
Although the sign of J\ is very sensitive to the O-D dis¬ 
tance in the crystallographic structure, a ferromagnetic 
Ji is only achieved in the DFT calculations of haydeeite 
for an O-D distance 20% shorter than the experimen¬ 
tally determined value of 0.98 A [20]. The origin of this 
discrepancy is not known. 

Diffuse inelastic magnetic scattering is clearly observed 
at temperatures above Xc, as shown in Fig. [3jb) at T = 
8 K. In contrast to the kagome antiferromagnets kapell¬ 
asite (S = 1/2) [21] and deuteronium jarosite (S = 5/2) 
[ 25 ], this scattering is clearly dispersive and mimics the 
spin-wave dispersion, thus resembling so-called param- 
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FIG. 4. (Color online) Cuts along Q of the magnetic scatter¬ 
ing x // (Q?X) at constant energies E as specified in the ferro- 
magnetically ordered phase of haydeeite o-MgCu 3 (OD) 6 Cl 2 
(closed circles) and the powder-averaged spin-wave calcula¬ 
tion for Ji — — 38 and Jd — +11 K (solid line), (a)-(b) Exci¬ 
tation energies E =s 1.0T0.015 and 1.4T0.015 meV measured 
with Ei = 3.55 meV at T = 1.5 K. (c)-(d) Excitation energies 
E — 8 ± 1 and 9 d= 1 meV measured with Ei = 16.6 meV at 
T = 2.2 K. 


agnons. Notably, this diffuse scattering seems to coex¬ 
ist with the spin waves below Xc, which can be inferred 
from the anomalously large Q-width of the scattering [see 
Fig-i a)]. Complementary measurements using the pow¬ 
der diffractometer WISH at ISIS show that the diffuse 
scattering is centered around Q = 0 and corresponds to 
short-range ferromagnetic correlations with a correlation 
length of about 35 A at T = 6 K (see Fig. [H]). The in¬ 
tensity of the short-ranged ferromagnetic correlations is 
roughly constant between T = 3.5 and 10 K. On further 
cooling, this intensity decreases, concomitantly with the 
reduction in %X (see inset of Fig. [5]). The behavior of 
the rather strong diffuse dynamic magnetic scattering in 
haydeeite is thus very different from what is expected for 
frozen-in magnetic disorder and also from critical scat¬ 
tering, for which the correlation length should diverge at 
Tc . The presence of strong short-range correlations close 
to a quantum phase transition where ferromagnetic order 
disappears is intriguing. 

We will now discuss the role of the Dzyaloshinski- 
Moriya interaction in haydeeite. The DM interaction is 
of current strong interest in kagome ferromagnets [Mb 
133] and kagome antiferromagnets above the saturation 
field [34], since an out-of-plane DM vector opens up gaps 
at finite energies in the spin-wave spectrum. This leads 
to a topological structure of the magnon bands and a 
finite Berry curvature in momentum space, resulting in 
the magnon Hall effect [35j [36] . The DM interaction is 
intrinsic to the kagome lattice, since the midpoint of the 
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nearest-neighbor bond lacks inversion symmetry [see Fig. 
[lja)]. In haydeeite, just as in the jarosite mineral fam- 
iiy [3l|, the Cu0 6 octahedra of the crystal structure are 
tilted with respect to the kagome planes, which allows the 
DM interaction to have both in-plane, D xy , and out-of- 
plane, D z , components m- We consider the ferromag¬ 
netic case with the assumption that the magnitude of the 
DM interaction is smaller than |Ji|, and use the experi¬ 
mentally determined exchange integrals (see above). For 
D xy = 0, the magnetic structure is a collinear ferromag- 
net. A finite value of D xy leads to a k = 0 umbrella struc¬ 
ture with the magnetization perpendicular to the kagome 
plane and where the opening angle of the umbrella in¬ 
creases with increasing \D xy \. A finite value of D z cre¬ 
ates gaps at finite energies of the spin-wave spectrum 
between the different branches, provided the spins have 
a component perpendicular to the kagome plane. How¬ 
ever, our neutron data are not compatible with the T 3 
IR that describes an umbrella-type structure with mag¬ 
netization perpendicular to the kagome planes. Such an 
umbrella structure should have very low intensity for the 
(001) and (002) magnetic Bragg peaks, whereas the data 
show strong intensity there (see Fig. [2|. Also, no finite- 
energy gaps are observed in the spin-wave spectrum. In 
other words, there is no experimental evidence that the 
DM interaction plays a significant role in haydeeite. 

We conclude that haydeeite appears to be a rare exam¬ 
ple of an isotropic S = 1/2 Heisenberg kagome ferromag- 
net. The ordered moment is strongly reduced to <0.2 /ib 
and coexists with diffuse dynamic short-range ferromag¬ 
netic correlations below Tc . We speculate that the re¬ 
duction of the ordered moment is due to quantum fluc- 



FIG. 5. (Color online) Energy-integrated diffuse scattering 
from haydeeite for various temperatures after subtraction of 
the low-T data. Inset: Magnetic susceptibility data \T vs. 
T of a-MgCu 3 (OD) 6 Cl 2 at H = 5 T showing ferromagnetic 
ordering at Tc — 4.2 K and a weak field cooling/zero-field 
cooling (FC-ZFC) splitting at lower temperatures. 


tuations in the proximity to the quantum phase transi¬ 
tion between the Heisenberg kagome ferromagnet and the 
Heisenberg cuboc2 antiferromagnet. Within this picture, 
a small increase in the third-neighbor across-hexagon an¬ 
tiferromagnetic exchange interaction Jd would push hay¬ 
deeite into the magnetically ordered cuboc2 phase, and a 
further increase of Jd or a slightly ferromagnetic second- 
neighbor interaction J<± would lead to a gapless chiral spin 
liquid state with spinon Fermi surfaces (CSL-A), which 
breaks reflection and time-reversal symmetries [181 . Such 
a state may have been observed in kapellasite m, where 
the Jd exchange is mediated by Zn 2+ rather than Mg 2+ 
ions. Studies of these compounds under pressure or with 
some partial Zn-Mg substitution may allow experimen¬ 
tal access to quantum spin-1/2 kagome systems at the 
border of a ferromagnetic instability. 
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